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二 阶 Dirichlet 问题 的 六 面体 等 参 元 逼近 
及 其 数值 积分 的 影响 * 


张 斐然 
(商丘 师范 学 院 数 学 系 ， 河 南 商丘 476000) 


摘 要 : 本 文 研究 的 是 三 维 曲 边 区 域 上 二 阶 椭圆 型 问题 的 等 参 有 限 元 逼近 。 利 用 等 参 变换 技巧 和 有 限 元 插 
值 理论 ， 本 文 证 明了 等 参 有 限 元 对 于 曲 边 区 域 边 值 问题 具有 最 优 收 和 敛 性 ， 保 证 了 等 参 有 限 元 数值 
格式 在 工程 应 用 中 的 可 靠 性 和 有 效 性 。 本 文 第 二 部 分 致力 于 构造 简单 有 效 的 数值 积分 公式 。 文 中 
给 出 的 四 点 、 五 点 、 六 点 和 八 点 六 面体 数值 积分 公式 ， 可 以 简化 单元 刚度 和 矩阵 和 荷载 向 量 的 计 
算 ， 并 严格 证 明了 利用 数值 积分 公式 的 等 参 元 数值 格式 仍然 保持 最 优 收 敛 阶 。 

关键 词 : 六 面体 等 参 元 ， 等 参 变 换 ， 数 值 积分 

分 类 号 : AMS(2000) 65N30; 65N12 中 图 分 类 号 : O241.82 文献 标识 码 : A 


有 限 元 方法 自 出 现 以 来 ， 由 于 其 强大 的 功能 和 相对 严格 完善 的 理论 ， 在 数学 界 和 工程 界 的 
各 领域 得 到 了 广泛 的 应 用 。 但 有 限 元 方法 处 理 曲 边 区 域 问题 时 ， 会 遇 到 网 格 和 区 域 不 完全 匹配 
的 问题 。 常 用 的 方法 有 两 种 :一 是 用 多 角形 区 域 来 逼近 ， 通 过 边界 处 网 格 的 充分 加 密 来 减 小 误 
差 4 。 但 这 要 面临 局 部 网 格 加 密 的 问题 ， 否 则 会 带 来 计算 量 的 巨大 增加 。 二 是 采用 等 参数 有 
限 元 方法 ， 这 种 方法 的 优点 是 可 以 保持 直 边 有 限 元 的 计算 量 ， 而 且 可 以 更 高 精度 的 逼近 区 域 的 
边界 。 但 由 于 引入 了 非 线性 变换 ， 给 有 限 元 格式 的 收敛 性 分 析 和 误差 估计 带 来 很 大 难度 。 另 
外 ， 由 于 变换 的 复杂 性 ， 程 序 的 实现 过 程 中 各 种 积分 的 计算 变 得 更 加 复杂 。 由 于 等 参 元 的 上 述 
优点 ， 使 得 该 方法 在 实际 中 (特别 是 工程 中 ) 得 到 很 多 应 用 ， 其 数学 文献 多 见于 在 二 维 问题 中 的 
应 用 区 9 引 ， 关 于 三 维 等 参 元 方法 的 分 析 和 计算 在 文献 中 并 不 多 见 。 

另外 ， 在 实际 计算 中 ， 单 元 刚度 和 矩阵 和 荷载 向 量 的 快速 计算 对 大 型 工程 问题 也 是 至 关 重 要 
的 一 个 环节 。 寻 找 单元 上 简单 易 算 的 数值 积分 公式 ， 成 为 一 个 有 意思 的 研究 课题 。 明 平 兵 等 在 
文献 [10] 中 研究 了 数值 积分 格式 对 旋转 Q 元 的 影响 ， 并 构造 了 简单 的 积分 公式 。 文 献 [11] 研 
究 了 数值 积分 公式 对 于 四 面体 等 参 元 的 影响 。 本 文 构造 了 应 用 于 六 面体 等 参 元 计算 的 四 种 简单 
有 效 的 数值 积分 格式 。 并 且 这 四 种 积分 格式 分 别 用 八 个 ， 六 个 ， 五 个 和 四 个 积分 点 ， 我 们 证 
明 ， 离 散 问 题 经 过 这 四 种 数值 积分 公式 的 近似 计算 后 ， 仍 然 收敛 于 连续 问题 的 解 ， 并 且 保 持 收 
敛 阶 不 变 。 因 此 ， 相 比较 而 言 ， 在 保持 整体 离散 精度 的 前 提 下 ， 积 分 公式 的 积分 点 越 少 ， 计 算 
就 越 简 单 。 
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2 ”六 面体 网 格 剖 分 和 等 参 有 限 元 


Oc R3 为 有 界 Lipschitz 区 域 ，T7 为 8 的 六 面体 剖 分 ， 其 所 有 单元 的 顶点 均 在 介 内 。 下 
面 我 们 通过 介绍 等 参 变换 来 引入 Q 的 六 面体 网 格 训 分 。Qk( 开 ) 是 天 上 对 每 一 个 分 量 的 指数 都 
不 超过 的 多 项 式 空间 ， 维 数 是 (k 十 1)"， 这 里 nn 是 空间 维 数 。Qk(K) 上 半 模 定义 为 


eons GL > 名 | az g > | Dmv(z)(e™ ) Paz) ue 


[v] = max sup || D™v(z)(e;™) || = max sup ge) | 
mR m Ee eS, : 1<i<n pek | DOzimm l 


其 中 e;™ = (eiei , 6i)。 
令 下 为 中 心 在 原点 ， 边 平行 于 坐标 轴 ， 边 长 为 2 的 正方 体 的 参考 单元 , (l < < 
8) 为 下 的 八 个 顶点 ， 其 中 
a1 = (-1,-1,-1), G2 = (1,-1,-1), a3 = (-1,1,-1), G4 = (1,1,-1), 
as = (-1, =1, 1), ag (1, —1,1), a7 = (-1, 1, 1), ag = (1, 1, 1), 


对 任 一 单元 K € TO， 定义 等 参 变换 Fk € (QÉ) A 
z= Fx(é = Daxhile md) viek, 


其 中 
PilE n, C) = (EEEL + mn + Gd) 
8 ? 


(Eom Gi) 为 正六 面体 的 第 ;个 顶点 的 自然 坐标 ， 构 成 Qi( 天 ) 的 一 组 基 。 容 易 验证 ，Fx 满 
Æ Fe (Gi) = a,x, 1<1<8. 
因此 ， 任 意 六 面体 单元 天 上 的 有 限 元 三 元 组 (K, Pe, UK) 定义 为 


Pr =spanf pi := pio Fr!, 1Si <8}, Ux ={aix,1<i<8}, 


所 有 这 样 的 单元 集合 (K, Px, Dx) 称 为 (1) 型 矩形 (K, Qi (K), Dg) 的 等 参 族 。 
KIR O 的 六 面体 单元 前 分 胞 为 


Th = {K : K = Fx(K), VK ETP}. 


3 ZB Dirichlet ABHA MAES TEE 
设 Q 是 Re" 中 有 界 开 子 集 ， 具 有 曲面 边界 9692， 考虑 二 阶 齐 次 Dirichlet 问题 
= 了 (9)， 
a(u,v) = fo = a430;u0,vdz, (1) 


f(v) = Jo fvdz, 
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其 中 ai e Ze(O)， je L2(Q)， 在 人 上 处 处 有 定义 ， 且 满足 椭圆 性 条 件 : 3 8 > 0， 使 得 
》 ayl > BD, 6&1l<i<n, vred. 
i,j=1 i=1 
我 们 用 (1) 型 六 面体 等 参 元 逼近 问题 (1)， 其 离散 空间 为 
Xp, = {un : Uh |KE Px } C C°(On), 


其 中 人 Qh = Uken K 一 般 来 说 £0, H Xr C C(O»); 因而 Xa C Hi(Q,)。 
令 


Vi = { vn E Xp: vn(b) = 0, b È 00 上 的 节点 }. 


显然 ， Vh C H! (Qh). 
AX OO, 和 59 非常 接近 ， 所 以 我 们 可 假定 ， 存 在 一 个 有 界 开 和 集合 使 ， 9 C 8， A, CO, X 
任意 的 h。 于 是 式 (1) 的 等 参 有 限 元 逼近 为 


| Kip, € Vi, 4 


人 z (2) 
Jo, 2o Gig Oitindjrnde = Jon funda, y Uh E Va; 
ij= 


其 中 Ga， 分别 是 ah f 到 锡 上 的 某 种 扩张 。 
现在 ， 在 等 参 元 逼近 问题 (2) 中 ， 以 数值 积分 代替 精确 积分 
| Ku, E Vn, 1E 


an(un, Vh) = (f, Un)h, V Un E Va, 


(3) 


其 中 


L 


n(Un, Up) = Det Denn ( oom vn) (b,x) }, 


K 


(F, un) = X. { yo K(fvn)(bi,K) } 


K 


考察 上 述 问 题 的 确定 性 ， 只 要 考察 数值 积分 节点 以 K BAAS KNOL. 
定义 1 车 存 在 与 无 关 的 常数 & > 0， 使 对 任意 的 v E Vh A 


dllvnll?,o, < on(vh, va), 


则 称 ap(-,-) 是 一 致 VW 椭圆 的 。 
4 SCHR [4] P §11 定理 1， 易 证 下 列 引 理 。 
引 理 1 车 an(:) 是 一 致 W, 椭圆 的 ， 则 存在 常数 C 与 九 无 关 ， 使 得 ;对 任意 的 立 e HQ), 


(Up, Wh) — Qh(Vh, W 
E E- c(i imf Ni valn + sup lã (vh, wa) — an (Vh, wa )| 
WhEVh walli, an 


t mp Bau Aatu], 
wh EVh llwallaon 
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qk 
a 


516 I # KR 学 


其 中 心 是 (3) 的 解 。 
我 们 用 本 kv 表示 函数 v 的 Pr- 插值，Hav 表示 函数 ,的 WW- 插值 ， 且 


Hxgv(ai kg) =v), Havlk =Ikv, Il6(6:) = 0(a). 


定理 1 设 {VW} 是 一 个 由 正则 (1) 型 六 面体 等 参 元 构成 的 有 限 元 空间 族 ，m < 3， 则 存 
HEC HAE, 对 任意 的 ve HQ) 


le — Havllmo, < Ch? ollz, m=0,1. 


证 明 因 n < 3, 2 > 3 > 8 HR) > COR) KENG) 型 六 面体 等 参 元 插值 误差 估计 


2 2 


式 
1_1,92. 
jv—Irgvlmark < O(h) a rh] olap K, m=0,1. 
此 处 p 二 4 二 2， 则 有 
lv 一 Txvlokg < < Ch2 «vie, KS Ch2 x lloll2, K, 


lv — Ixvh,x < Chrlvlz x < Chx|lvllo,x, 


AA Tune = Uku H — 0K c ,成立 ， 故 我 们 有 
lv- Tolo, = (Yo le- Telk)? <cn( E lol)” = Chole 
KETh KETh 


定理 证 毕 。 


4 ”数值 积分 的 影响 

本 节 将 考察 数值 积分 对 六 面体 等 参 有 限 元 逼近 问题 (1) 的 影响 ， 考 虑 四 种 数值 积分 格式 ， 它 
们 在 参考 单元 上 分 别 用 八 个 、 六 个 、 五 个 和 四 个 积分 点 ， 对 这 四 种 数值 积分 格式 ， 我 们 将 证 明 
EM H) 误差 估计 。 

首先 ， 定 义 逼 近 问 题 。 

问题 1 Kun E Vr, Falun, v) =(f,v), Yv E Vro 

在 参考 单元 上 定义 数值 积分 格式 


I plêd($) ~ Sap, wi>0, hek, 1<i<I. 
数值 积分 误差 泛 函 记 为 


Ex(e)= f iE ened, 


i=l 


Be) = f, padt- SOP), Erle) = BelGe) 


i=l 


AQ = {1,&,m,C,&n,EC,nC,EnC}。 假 定数 值 积 分 格式 对 QT, BMF RUA: 
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格式 1: I=8 
{i}21 = 1,1,1,1,1,1,1,1, 
{i} = (-1,-1,-1), (1, —1, —1), (-1, 1, -1), (1, 1, —1), 
(—1,-1, 1), (1, —1, 1), (—1, 1,1), (1,1,1), 
格式 2: I=6 
{2:}& = 1,1,1,1,2,2, 
{6:}$.1 = (-1,-1,-1), (1, —1, —1), (-1,1,-1), (1,1, -1), (0, —1, 1), (0,1, 1), 
格式 3: I=5 
{@,}8_, = 1,1,2,2,2, 
{6}. = (~1, —1, —1), (1, —1, —1), (0, 1, —1), (0, —1, 1), (0, 1, 1), 
格式 4 I=4 
{2i} = 2,2,2,2,  {8;}4., = (0, —1, —1), (0,1,~1), (0, —1, 1), (0, 1, 1), 
{2i} = 2,2,2,2,  {6;}4., = (-1,0, -1), (1,0, —1), (-1,0, 1), (1,0,1), 
{Gi}in1 = 2,2,2,2,  {6;}4., = (-1, —1,0), (1, -1,0), (—1, 1,0), (1, 1,0), 
注 格式 2 与 格式 3 仅 考 虑 了 其 中 一 种 情形 ， 其 它 几 种 情形 可 对 称 得 到 。 
现在 ， 把 数值 积分 格式 应 用 到 有 限 元 方程 


3 
a,(u,v) = >> Qx( > oujBiu8iv ) 及 (fv)h= J Qr(fvn), 


KETk i,j=l1 KETk 


则 问题 1 转化 为 求解 如 下 问题 。 
问题 2 Kur E Vn, 使 得 arlun, v) = (fv)h, Vu € Vino 
引 理 2 修正 变 分 形式 on(.,:) 是 VV 椭圆 的 ， 即 存在 与 无关 的 常数 CC， 使 得 


clvnlli,g, < an(vh, va), Von € Va. (4) 


证 明 ”由 文献 向 中 811 定 理 3 知 ， 对 积分 格式 1 至 格式 4 来 说 ，(4) 式 显然 成 立 。 
Oi SCHR [11] 中 引 理 5.2 和 引 理 5.3 的 证 明 ， 可 得 如 下 结论 。 
引 理 3 设 已 给 定 一 个 (1) 型 矩形 的 正则 等 参 族 (天 )， 参 考 元 忘 上 数值 积分 公式 满足 


E(g)=0, VGEQi(K), 


则 存在 常数 CC 与 单元 KK 无 关 ， 使 得 对 任意 的 aij € Wl%(K), we Pes #8 


3 3 
| Z Bx (ai50;v0,w) | < Chic ( Y lolio, ) lolar holsr. 
1 i, j=1 


ij= 
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引 理 4， 设 已 给 定 一 个 (1) 型 矩 形 的 正则 等 参 族 (入 )， 参 考 元 并 上 数值 积分 公式 满足 
E(g)=0, Vee Qi (K), 
则 存在 常数 C ATK 无关， 使 得 对 任意 的 fe W29(K), ve Pe, g > 2， 有 


1 


|Ex(fv)| < Chg (mean(K))?~* |i flia, lloll,- 
定理 2 K{V,} Æ (1) 型 矩形 的 等 参 族 ， 若 满足 
(i) EWER: FEBR o My PERN K, KEHNA 


h 
<a, YK <7 <1; 
hg 


(区 ， 忘 上 数值 积分 公式 满足 ， 对 任意 的 六 Q1(K), E(G) = 0; 
(iii) 设 癌 是 一 个 开 集 ， 满 足 : 2900, n CH, HERR NA 是 (1) 型 六 面体 构成 
的 Q 近 似 区 域 ， 且 ww 和 oij 存在 一 个 扩张 入 和 Gi;， 定义 在 全 上 满足 
ŭe H?(Ñ), Gi EWHOO), 1<ij<3 
E 3 
f=- DD 8; (igi) E W™(Ñ),  g = 2, 


i,j=l 
则 存在 常数 C， 使 得 
3 
lü — urlio, < Ch( lalla + 2 lastly allze + fl ). 
i,j=l 
证 明 ”我 们 应 用 引 理 1， 就 有 


lalat, wa) 一 ar (Ini, wa )| 


lë — unlon, < O(a Taüll + sup > 
a ate ja wn) — fr(wa)| ), 
wnreVn lwallso, 
Hh ma BRR OVA, AA < 2， 所 以 H29) 6 OÀ), MERA MN. 
由 定理 1， 我 们 有 


其 次 ， 利 用 引 理 2， 引 理 3 K Schwarz FER, RNA 


|| Trill,o, < chilly 5 


d(T, wa) — an (Mnt, wn) 


3 
= 5 5 Ex (GijOiIn G0; wn), (因为 Qij (bik) T Gij (bi,K)) 


KET, 1,j=1 


ch( > lësh oa )( 5 Mälk) onlay 


i, j=1 KET, 


IA 


3 
Ch Yo az lllz allwn. 


i,j=1 


[人 
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|@(IIpt, wr) 一 an (TInt, wa )| Be ss $ 
< Ch > laszl ,0,all@lle,9- 


ij=l 


sup 
wn EV walli, an 


由 于 数值 积分 节点 bi,kx EQ, PRU f(bi,k) = f(bi,k)。 这 样 ， 我 们 利用 引 理 4 得 


a(t, Wr) 一 fh (wh x 
sup [ats tn) = fawn)! _ SO sup Ex(fun) 
wh EVh walli, a KET, WhEVh 


Ch( >》, mean(K))* “filia = Chl|fllya0,. 


KET, 


IA 


定理 证 毕 。 
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The Approximation by Hexahedral Isoparametric Element and the 
Effect of Numerical Integration of Second Order Dirichlet Problem 


ZHANG Fei-ran 


(Department of Mathematics, Shangqiu Teachers Colloege, Shangqiu Henan 476000) 


Abstract: In the context of numerical simulations, three-dimensional domains with curved boundaries 
are difficult to deal with, since many numerical schemes can not capture the boundary information ef- 
ficiently. In this paper, we study the hexahedral isoparametric element method for elliptic problems on 
domains with curved boundaries. The optimal convergence of the discrete problem to the continuous 
problem is obtained upon using the isoparametric transformation and scaling techniques. To simplify 
the computation of elementary matrices and load vectors, we develop four numerical integration for- 
mulas on hexahedrons. The optimal convergence of the isoparametric element method using our four 
integration formulas are also obtained. 
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